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REPORT  SUMMARY 

At  the  conclusion  of  this  quarter  the  following  milestones 
were  agreed  upon  by  the  University  of  Illinois  and  the  Burroughs 
Corporation: 

Event  Date 


Debugged  Prototype  PE  Oct.  '69 

Debugged  Production  PE  Feb.  '70 

CU  Ready  for  Integration, 

Start  System  Integration  Feb.  '70 

System  Confidence  Demonstration        Apr.  '70 

Ship  One-Quadrant  System  Aug.  '70 

Final  Acceptance  Test 

at  University  of  Illinois  Nov.  '70 

The  semi-conductor  memories  to  be  produced  by  Fairchild  are 
two  weeks  to  one  month  behind  schedule.   The  memory  chips  for  the 
prototype  processing  element  memory  appear  to  be  sound  and  operable, 
however,  Fairchild  is  experiencing  yield  problems.  Fairchild  is 
currently  developing  a  new  set  of  procedures  and  equipment  to  test 
memory  chips  at  the  wafer  level. 

Texas  Instruments  is  currently  meeting  Burroughs'  delivery 
schedule  in  providing  the  dual-in-line  packaged  ECL  circuits  for 
ILLIAC  IV. 

The  PDP-9  computer,  which  will  be  attached  to  the  PE  exerciser 
to  provide  the  ability  to  automatically  diagnose  errors  in  the  PE  and 
CU  boards,  was  delivered  to  the  University  in  August.   Components  for 
the  interface  to  the  PE  exerciser  have  arrived  at  the  University  and  the 
interface  is  currently  under  construction.   Software  development  for  the 
FDP-9  is  proceeding  on  schedule. 

A  new  version  of  the  "ILLIAC  IV  Systems  Characteristics  and 
Programming  Manual"  has  been  delivered  by  Burroughs  to  the  University 
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and  distributed  to  the  appropriate  personnel.   The  University  and 
Burroughs  are  working  together  to  gather  comments,  criticisms  and 
corrections  to  he  incorporated  into  the  next  update  of  the  manual  to  he 
published  in  January ,  1970- 

FOURUM,  the  ILLIAC  IV  user's  group,  continues  to  grow  and  the 
present  membership  is  at  l4o.   The  distribution  for  other  research 
documents  is  in  the  process  of  being  automated  and  all  mailing  lists  are 
being  centralized  into  one  file. 

Burroughs  has  agreed  to  accept  an  incremental  conversion  of 
the  present  contract  to  a  fixed  price  contract.   Target  dates  for 
specific  tasks  which  are  to  be  converted  to  fixed  price  are  currently 
\ander  negotiation. 

The  project  initiated  last  quarter  to  permit  Burroughs'  ALGOL 
to  be  compiled  and  executed  on  a  Control  Data  Corporation  60OO  series 
machine  has  been  partially  completed.   A  preprocessing  program  that 
translates  Burroughs  ALGOL  and  XALGOL  into  CDC  ALGOL  is  currently 
operational.   A  more  comprehensive  program  should  be  available  next 
quarter. 

The  building  plans  for  the  Center  for  Advanced  Computation 
were  approved  by  the  President  of  the  University  of  Illinois  and  the 
Board  of  Trustees.   Bids  for  general  construction  have  been  let. 
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1.   HARDWARE 

1.1  Logic  Simulation  and  Diagnostics 

1.1.1  Logic  Simulation 

1.1.1.1  PE  Simulator 

The  CU  Card  Simulator  Generator  System  has  been  extended  and 
applied  to  the  generation  of  EE  Card  Simulators.   The  net  lists  (wire 
lists)  of  35  PE  card  types  vere  compiled  from  Burroughs'  logic  diagrams. 
Error  messages  from  the  programs  of  the  Generator  System  helped  in 
removing  some  errors  in  translation  and  in  card  punching  of  the  net 
list  as  well  as  in  the  logic  diagrams.   The  connector  pin  assignment 
in  the  diagrams  was  compared  with  Burroughs'  assignment  list  which  was 
being  used  in  the  generation  of  the  backplane  wire  list.   All  of  the 
HI  card  net  lists  have  been  released  for  simulation  and  test  generation. 

1.1.1.2  CU  Simulator  System 

The  CU  Card  Simulator  Generator  System  is  now  complete ^  with 
the  facilities  for  automatic  sequential  execution  and  history  generation. 

Some  experiments  of  CU  Card  Simulation  were  made  and  the 
system  was  verified. 

The  system  will  be  used  for  debugging  the  logic  of  CU  cards 
and  for  generating  Bad  Card  Simulators. 

The  design  of  the  CU  Section  Simulator  Generator  System  was 
resumed  in  September.   The  System  will  be  implemented  in  the  next 
quarter  and  used  for  checking  the  design  of  the  Control  Unit. 

1.1.2  Card  Test  Generation 

1.1.2.1  Card  Test  Generator 

In  this  period,  a  final  version  of  the  Card  Test  Generation 
System  has  been  completed  and  debugged. 
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To  allow  more  flexibility  and  to  meet  additional  requirements 
(e.gv  initialization  of  the  latches  present  in  some  "boards  and  holding 
some  signal  values  at  a  constant  during  the  simulation  instead  of 
varying  them  in  a  random  manner),  a  new  program  has  been  incorporated 
into  the  system. 

It  is  expected  that  the  capabilities  of  the  system  will  be 
extended  in  the  near  future  to  facilitate  failure  isolation  in  addition 
to  failure  detection. 

1.1.2.2  Card  Test  Translation 

The  two  main  functions  of  the  Card  Test  Translator  are:   (l)  to 
print  out  the  test  patterns  and  the  list  of  failures  associated  with 
patterns  in  an  adequate  form,  and  (2)  to  translate  the  patterns  and 
failure  information  into  PEXTAP  (EE  Exerciser  Test  Assembler)  language. 

The  first  function  has  been  completed  and  the  second  function 
is  being  debugged.   The  program  will  be  completed  in  the  next  quarter 
with  some  extension  of  the  functions. 

1.1.2.3  Generation  and  Evaluation  of  Card  Tests 

More  than  half  of  the  35  PE  cards  have  been  simulated  for  test 
generation.   In  most  cases,  a  set  of  patterns  which  can  detect  all 
failures  in  the  DIL  packages  on  the  card,  was  obtained  in  a  short  running 
time  (between  10  seconds  and  2  minutes)  on  the  B5500.   The  rest  of  the 
time  required  on  the  computer  was  about  10  minutes  including  simulator 
compilation  and  random  pattern  generation. 

Also,  a  few  CU  cards  were  simulated  for  test  generation  during 
this  period.   All  failures  were  detected  on  one  of  the  cards.   The  input 
partitioning  and  the  initialization  of  the  storage  states  seem  to  be 
very  important  for  efficient  test  generation. 

Burroughs  has  compiled  a  list  of  possible  failure  modes  in 
DIL's  by  analysis  of  the  circuits.   Work  is  continuing  at  the  University 
to  correlate  their  results  to  the  failure  modes  assumed  in  the  Card  Test 
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Generator.  It  appears  that  practically  all  failure  modes  defined  by 
Burroughs  can  be  covered  as  single  or  double  failures  in  the  Card  Test 
Generator. 

1.1.3  PE  Diagnostics 

The  path  test  generator  system  was  used  to  generate  a  revised 
test  program  for  the  new  prototype  PE.   Combinational  tests  have  not  yet 
been  revised  for  the  new  design. 

One  of  the  last  programs  of  the  Control  Logic  Test  Generator 
System  is  being  coded.   This  program  will  detect  and  solve  conflicts  in 
the  assignment  of  values  to  input  signals. 

1.1.4  PEX  Computer  and  Supervisor  System 

Delivery  of  a  Digital  Equipment  Corporation  PDP-9/L  computer 
to  be  used  as  a  PEX  controller  was  made  at  the  end  of  August.  Charac- 
teristics of  this  system  are: 

1.5  (J.sec  cycle  time 

18  bit  word  length 

16k  words  of  core  storage 

Hardware  multiply  and  divide 

3  DECtape  transports 

1  IBM  compatible  magnetic  tape  transport 

Coding  and  assembly  of  the  previously  defined  PEX  supervisor 
is  progressing  satisfactorily.   Early  hands-on  experience  with  the 
computer  was  gained  through  coding  and  assembling  of  several  utility 
routines  to  speed  development  of  the  PEX  supervisor.   Typical  of  these 
is  a  magnetic  tape  handler  which  permits  listings  to  be  written  on  a 
magnetic  tape  whose  contents  are  then  printed  using  the  high  speed  line 
printer  of  the  Burroughs  B5500. 

The  interface  to  the  PEX  will  be  made  by  modifying  the  second 
paper  tape  reader  channel  on  the  PEX,  and  then  constructing  an  interface 
to  the  PDP-9  which  will  accept  data  words  from  the  PDP-9's  data  channel 
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and  serialize  them  into  characters.   Design  of  the  interface  was 
accomplished  during  August.  A   misunderstanding  "between  Burroughs'  and 
the  University's  logic  designers  ahout  the  exact  details  of  the  inter- 
face required  that  small  modifications  he  made  in  mid- September. 

Actual  construction  of  the  interface  has  been  delayed  awaiting 
arrival  of  the  necessary  components.   It  is  expected  that  the  interface 
construction  will  be  completed  at  approximately  the  same  time  as  the 
programming^  and  that  the  entire  system  will  be  delivered  to  Burroughs 
in  time  to  check  out  the  first  production  PE's. 

1.2  Design  Automation 

1.2.1  Bost  Processor  Program 

The  Post  Processor  (Printed  Circuit  Wiring  Analysis)  Program 
is  being  used  for  production  of  the  CU  printed  circuit  boards.   The 
program's  output  has  been  useful  as  an  additional  tool  for  updating  and 
completing  the  PC  Artwork.   The  need  for  substantial  manpower  to  manually 
check  the  automatically  produced  artwork  has  been  eliminated. 

1.2.2  CU  Printed  Circuit  Layout 

The  University  and  Burroughs  have  completed  arrangements  for 
doing  the  DA  for  the  CU  PC  layout.   Approximately  12  hours  of  processing 
is  done  on  the  University's  B5500  system  per  PC  board.   The  University's 
system  can  provide  better  turnaround  and  throughput  for  this  larger 
amount  of  processing  and  save  additional  billing  to  the  project.   It  is 
planned  to  complete  this  effort  during  the  next  quarter. 
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2.   SOFTWARE 


2.1  Operating  System 

Several  modifications  were  made  to  the  operating  system  to 
provide  for  the  first  pass  October  system.   A  second  job  parser  was 
designed  and  is  being  debugged.   This  job  parser  is  a  simpler  module 
than  the  first  and  should  be  more  reliable  during  early  debugging  of  the 
first  system.   The  loader  functions  have  been  subsumed  by  the  program 
collector  and  execution  monitor.   Since  multiple  relocation  and  overlay 
loading  would  not  be  required  in  the  first  system,  the  program  collector 
was  modified  to  produce  a  relocated  program  segment  suitable  for  direct 
loading  into  the  ILLIAC  IV  memory.   Thus,  the  loader  was  made  unneces- 
sary and  the  read  operation  needed  to  place  the  program  in  the  array  was 
assigned  to  the  execution  monitor. 

The  hardware  supervisor  intrinsics  were  frozen  and  the  imple- 
mentation of  buffer  files  on  the  B65OO  has  begun.  Contrary  to  expecta- 
tions. Burroughs  will  not  include  the  buffer  file  facility  in  its  first 
Master  Control  Program  (MCP) . 

The  execution  monitor,  data  processor,  and  job  partner  depend 
heavily  on  the  MCP.   Due  to  the  moving  target  nature  of  the  MCP,  the 
specification  of  critical  sections  of  each  of  these  modules  was  not  made 
until  mid -September.   It  was  decided  to  freeze  the  design  on  the  basis 
of  current  knowledge  of  the  MCP  and  to  modify  the  delivered  MCP  to 
represent  the  situation  on  which  the  design  was  frozen.   It  is  expected 
that  this  may  entail  an  extensive  amount  of  debugging  time  when  the 
B650O  is  first  delivered. 

Parts  of  the  system  have  been  coded  in  SIMULA  so  that  their 

activity  can  be  simulated  and  debugged  before  the  arrival  of  the  B65OO. 

Currently,  simulation  activity  is  concerned  with  disk  allocation  and 
quadrant  execution. 
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2.2  Translator  Writing  System 

Most  of  the  steps  of  the  system  described  last  time  have  been 
implemented  and  effectively  debugged.   BIIF2FPL/TWS,  FPL2PAR/TWS, 
PAE2ALG/WS,  and  the  parser  files  are  all  functional.   TWIMCLE/DISK  is 
working,  except  for  a  few  of  the  TVTIWKLE  constructs;  work  on  ISL/DISK 
is  just  beginning.   Use  of  the  new  system  has  already  begun.  From  the 
results  obtained  so  far,  the  expectations  for  the  new  system  have  been 
met.   The  interpretive  parser  is  very  similar  in  all  respects  to  the 
previous  one,  but  the  executable  parser  is  significantly  smaller.   The 
executable  parser  is  small  enough  (13K  core  estimate  for  executable 
TESLA  compiler  compared  with  lOK  for  the  interpretive  version  and  18k 
or  20K  for  the  old  executable  version)  that  it  has  become  a  very  desir- 
able choice  for  small  to  medium-sized  languages.  A   further  advantage 
is  that  three  different  pieces,  the  filling  of  the  parsing  tables,  the 
lookahead  procedure,  and  the  Floyd  production  procedure  calls,  can  be 
made  executable  or  interpretive  independently. 

2. 3  Compilers  and  Translators 

2.3.1  TRANQUIL 

Progress  on  the  TRANQUIL  compiler  was  slowed  during  this 
quarter  due  to  personnel  changes  and  vacations.   Brisk  activity  on  the 
compiler  was  resumed  during  the  last  two  weeks  of  the  quarter.  While 
coding  and  debugging  of  the  TRANQUIL  basic  subset  due  for  completion  in 
October  continues,  development  and  design  effort  has  been  centered  on 
implementation  of  explicit  l/o  statements  for  the  language. 

The  implementation  of  procedures  in  TRANQUIL  has  been  com- 
pleted on  time.   Some  slippages  are  apparent  in  the  coding  and  debugging 
of  the  basic  subset.   It  was  hoped  to  have  coding  completed  by  the  end 
of  October,  but  November  is  now  the  target  date.   Debugging  will  take 
longer  and  will  be  a  function,  in  part,  of  the  user's  perseverance  in 
utilizing  what  is  available. 
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2.3.2  GLYPNIR 

Work  on  GLYPNIR  proceeded  slowly  during  this  quarter  because 
of  a  reduction  in  staff.  Even  so,  several  additions  and  extensions  to 
the  compiler  were  made.  The  compiler  now  accepts  numeric  literals, 
including  floating  point  numbers  with  exponents  and  base  designators. 
Also,  partial  word  operators  (similar  to  those  of  BURROUGHS  EXTENDED 
ALGOL),  the  AS  construct  (which  allows  the  programmer  to  make  a  GLYHHR 
register  synonymous  with  an  ILLIAC  IV  hardware  register),  and  several 
new  control  card  options  have  been  added  which  will  make  debugging 
easier. 

Work  is  also  progressing  on  the  debugging  of  system  sub- 
routines such  as  sine,  cosine,  tangent,  natural  log,  and  square  root. 
These  system  subroutines  will  be  introduced  during  the  early  part  of 
the  next  quarter. 

2.k     Assembler 

Juily  and  August  were  spent  in  fixing  the  final  design  for  the 
macro  assembler.   Implementation  of  the  macro  assembler,  according  to 
this  design,  began  in  September. 

2.5  Utilities 

The  utilities  effort  commenced  during  this  quarter.   The 
initial  effort  is  concentrated  on  developing  two  formatted  data  process- 
ing programs.   One  will  take  data  from  external  sources  and  put  it  on 
the  B65OO  disk  in  the  form  required  by  the  operating  system  for  trans- 
mission to  the  ILLIAC  IV  disk.   The  other  program  performs  the  inverse 
transformation:   it  takes  a  file  which  v/as  transmitted  to  the  B65OO  disk 
from  the  ILLIAC  IV  disk  and  formats  the  data  for  external  representation. 

Other  areas  under  consideration  are : 

(a)  debugging  facilities 

(b)  graphic  display  software 

(c)  software  to  permit  user  measurement  of  program 
performance  with  an  eye  to  hardware  usage  and  time 
optimization. 
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2.6  B^^OO  Operation 

The  present  status  of  the  MCP  is  Mark  IX,  LEVEL  51.1^.   In 
J-uly  175 '16  hours,  in  August  I'jk.'JO   hours  and  in  September  202.39  hours 
of  processor  time  was  logged. 
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3.   APPLICATIONS 


3.1  Numerical  Analysis 

3.1.1  Monte  Carlo  Transport 

A  two-dimensional  Monte  Carlo  Transport  program  is  ciirrently 
being  written  in  GLYPNIR.   Straight  storage  is  used  for  the  variable 
storage  scheme.   The  program  will  accommodate  approximately  6k  x  6k 
spacial  cells,  and  only  I80  x  6k   particles,,  but  appropriate  buffers  are 
being  provided  so  that  particles  can  be  streamed  in  and  out  from  the 
disc.   The  coding  is  more  than  50  percent  complete. 

3.1.2  Mathematical  Subroutines 

The  following  subroutines  have  been  completed  and  simulated. 


Precision 

Function 

in  Octa 

15  ~ 

1  Digits 

sine 

16 

cosine 

15  ~ 

16 

tangent 

15  ~ 

16 

cotangent 

15  ~ 

16 

arccosine 

15  ~ 

16 

arc sine 

15  ~ 

16 

arctangent 

15  ~ 

16 

natural  logarithm 

Ik   ~ 

16 

exponential 

Ik   ~ 

16 

square  root 

15  ~ 

16 

The  subroutines  are  entered  with  6k   arguments  and  they  return 
6k   answers.   Computer  Approximations,  by  Hart  [1],  was  the  source  for  the 
approximating  functions.  Accuracy  was  determined  by  comparing  results 
from  the  ILLIAC  IV  simulator  with  double  precision  results  from  the 
360/75.   Corresponding  routines  for  32  bit  operands  will  be  provided  as 
they  are  needed. 
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3.1.3  Automatic  Solution  of  Initial  Value  Problems 

An  investigation  was  initiated  this  quarter  to  determine  the 
feasibility  of  constructing  a  general  purpose  program  to  solve  initial 
value  problems  of  the  form: 

f=irw  (1) 

where  A  may  be  a  non-linear  operator  involving  partial  derivatives  with 
respect  to  spacial  variables.   Especially  encouraging  in  this  investi- 
gation is  the  success  the  University  has  had  with  automatic  integration 
programs  to  solve  large  sets  of  ordinary  differential  equations  of  the 
form: 

f=?(?)  (2) 

at 

2  /  2 

It  can  be  shown  for  certain  operators  A  (for  example  S  /^x  ) 

that  many  of  the  successful  finite  difference  schemes  for  solving  (l) 
may  be  derived  by  first  discretizing  the  spacial  variable,  forming  a 
large  set  of  ordinary  differential  equations,  and  then  solving  the 
ordinary  differential  equations  using  a  conventional  scheme.   Consider- 
ing accuracy  and  stability  there  is  no  difference  in  the  approach.   This 
suggests  that  we  may  capitalize  on  the  success  of  the  automatic  integra- 
tion of  ordinary  differential  equations  in  the  integration  of  partial 
differential  equations. 

The  two  pressing  questions  which  must  be  investigated  are: 
(a)  Can  the  spacial  mesh  be  refined  independently  if  the  algorithm  for 
solving  the  ordinary  differential  equations  has  certain  properties? 
and  (b)  Will  it  be  possible  to  first  approximate  A  by  another  operator 
A*  which  will  insure  that  discontinuities  in  the  solution  (i.e.  shocks) 
will  not  completely  destroy  the  numerical  solution?   If  these  questions 
can  be  answered  affirmatively,  it  may  then  be  practical  to  develop  an 
automatic  solution  program. 
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3.1.4  Eigenvalue  Problems 

The  program,  for  Jacob! 's  method  for  eigenvalues  of  symmetric 
matrices  with  a  "Triangular  Skew"  storage  scheme,  has  been  written  and 
is  now  in  the  process  of  being  simulated  and  debugged. 

The  problem  most  apparent  at  the  present  time  is  finding  a 
suitable  method  of  keeping  track  of  the  elements  in  a  most  efficient 
way.   Since  in  each  transformation  of  the  matrix  only  those  row- 
elements  change  their  locations  whose  index  n  <  #  of  transformations, 
the  method  described  in  QPR  Document  #219  [2]  seems  inefficient,  since 
it  calculates  all  locations  of  all  elements. 

The  method  being  developed  sees  the  location  of  the  elements 
of  one  row  as  a  function  of  the  locations  of  the  elements  of  the 
previous  row  with  the  adjustment  that  each  row  has  one  element  less 
than  the  previous  one.   This  method  also  involves  a  minimum  of  element- 
manipulation. 

3.2  Linear  Programming 

Major  projects  for  the  Linear  Programming  (LP)  group  this 
quarter  have  been  the  SSK  simulation  of  LPS  (small  linear  programming 
system)  and  an  investigation  of  the  gradient  projection  method  as  a 
technique  for  LPL  (large  linear  programming  system). 

One  of  the  significant  benefits  that  linear  programming  by 
the  simplex  method  will  receive  from  ILLIAC  IV  is  in  the  realm  of 
"multiple  pricing".   Behind  this  concept  lies  the  fact  that  no  computer- 
implemented  LP  system  can  afford  to  examine  the  entirety  of  a  large 
constraint  matrix  (A-matrix)  each  iteration.   Typically,  reduced  prices 
(d^'s)  are  computed  for  only  one  portion  of  the  A-matrix  columns 
(partial  pricing).  A   smaller  number  (about  5)^  with  the  most  favorable 
prices,  are  selected  for  pivot  evaluations.   This  is  multiple  pricing. 
Updating  such  a  set  of  vectors  is  time-consuming,  so  suboptimization  is 
performed  with  this  set;  that  is,  as  many  of  these  vectors  as  possible 
are  introduced  into  the  basis  before  reconsidering  the  remaining 
A-matrix  columns. 
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The  massive  "crunching"  capability  of  the  ILLIAC  IV  allows  the 
entire  A-matrix  to  he  scanned  in  determining  the  pricing  vectors.   A 
major  iteration  consists  of  pricing  all  columns  of  the  A-matrix  and 
selecting  and  transforming  those  selected  for  multiple  pricing.   LPS 
handles  up  to  127  (vs.  1-10  on  a  conventional  machine)  of  these  vectors 
for  suhoptimization.   The  suboptimization  steps  are  called  minor 
iterations.   The  pivot  element  is  determined  for  each  of  these  trans- 
formed vectors;  the  one  producing  the  greatest  improvement  in  the 
objective  function  is  selected  for  pivoting.   When  the  remaining 
transformed  vectors  produce  insufficient  objective  function  improvement, 
a  new  major  iteration  is  initiated. 

The  pricing  and  candidate  selection  portion  of  the  major 

iteration  (DSCHQZ)  has  been  successfully  simulated  this  quarter. 

Consisting  of  three  sequentially  executed  steps,  DSCHUZ  (l)  calculates 

the  d-'s  for  the  entire  A-matrix  and  selects  an  optimum  sub set- -usually 
J 

512--of  these,  (2)  refines  this  set  to  N  <  127  vectors  and  (3)  executes 
an  inter-PE  sort  on  the  selected  values.   The  result  is  a  column- sorted 
array  containing  tags  and  cost  values  for  the  selected  vectors.   Then 
the  N  most  probable  vectors  to  enter  the  basis  can  be  column- sequentially 
processed  in  minor  iterations.   Simulation  tests  of  other  LPS  routines 
will  be  carried  out  as  machine  time  becomes  available. 

Investigation  of  gradient  projection  as  a  solution  technique 
for  LPL  (large  linear  programming  system)  has  produced  a  preliminary 
recursion  algorithm  for  the  method.   It  takes  a  form  which  is  similar 
in  concept  and  im.portance  to  the  recursion  introduced  by  the  product 
form  inverse  in  place  of  the  explicit  inverse  of  the  revised  simplex 
method.   Studies  of  the  vector  densities  produced  and  the  resulting 
storage  requirements  are  somewhat  alarming  at  present;  modifications 
axe   being  considered  as  the  investigation  continues. 

3. 3  Long  Codes 

During  this  quarter,  the  identification  of  a  stationary 
dynamical  system  with  plant  and  observation  noises  was  attempted.  A 
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least-squares  estimation  yielded  an  estimation  error  that  remains 
bounded,  but  it  did  not  provide  an  upper  bound  on  the  size  of  this 
error  in  the  limit. 

A  stochastic  approximation  algorithm  was  constructed  based 
on  some  theorems  by  Dvoretzky  [3]  and  by  Venter  [h].      It  was  shown 
that  the  proposed  algorithm  yields  an  estimation  error  that  approaches 
zero  in  the  mean- square  sense. 

3.^  Radar  Data  Processing 

During  this  period,  work  was  completed  on  the  documentation 
of  the  ILLIAC  IV  Kalman  Filter.  A   complete  write-up  has  been  printed 
and  is  available  as  ILLIAC  IV  Document  #222,  entitled  "A  Kalman-Filter 
Tracking  Program  for  ILLIAC  IV"  [5]. 

A  complete  subroutine  package  for  doing  Fourier  transforms  is 
being  written  for  IIIjIAC  IV.   It  uses  the  Cooley-Tukey  algorithm  and 
can  handle  data  sets  which  can  vary  from  8  to  ^096.   The  algorithm  is 
programmed  in  a  32-bit  floating  point  mode  for  complex  numbers,  where 
the  real  and  imaginary  value  of  each  reside  in  one  64-bit  PE.   This 
program  is  expected  to  be  running  by  the  next  quarter  and  complete 
documentation  will  be  started. 

3. 5  Seismic  Signal  Processing 

Several  programs,  previously  written,  were  modified  to  make 
them  more  flexible.   One  research  assistant,  working  on  these  programs, 
was  beginning  to  convert  these  B5500  ALGOL  programs  to  ILLIAC  IV 
assembly  language  xantil  the  draft  board  assigned  him  to  a  different 
area  of  responsibility. 

Research  into  multi-channel  filter  design  has  continued  during 
this  quarter. 

Since  the  data  received  in  seismic  recording  is  all  real, 
rather  than  complex,  a  modification  of  the  Fast  Fourier  Transform 
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program  to  accept  two  streams  of  real  data  has  started.  This  program 
has  "been  postponed  until  the  fourth  c[uarter  when  additional  personnel 
will  "be  acquired. 

3.6  Conversion  of  Burroughs  ALGOL  to  CDC  60OO  Series  ALGOL 

A  preprocessing  program  that  translates  Burroughs  ALGOL  and 
XALGOL  into  CDC  ALGOL  is  operational.   Both,  the  ILLIAC  IV  assembler 
and  simulator,  have  heen  properly  processed  by  this  program. 

Certain  Burroughs  constructs,  such  as  SCAN  and  REPLACE 
statements,  are  changed  into  function  calls.   All  references  to  input- 
output  are  also  changed  into  function  calls.   The  necessary  subroutines 
must  be  written  for  the  CDC  machine.   It  is  expected  that  work  should 
be  completely  finished  by  the  end  of  November  1969* 

Programs  will  run  20  percent  faster  on  the  CDC  machine  than 
previously  estimated. 

3.7  ILLIAC  IV  Education  and  Documentation 

3.7.1  cs  491-D 

The  graduate  course  in  Computer  Science   "Architecture, 
Applications  and  Languages  for  a  Parallel  Computer"  is  being  offered 
again  this  fall.   The  course  outline  is  as  follows: 

I.  Machine  Design  .  .  .  (Walt  Heimerdinger)  ....  6  classes 

Conventional  Machines 1 

Unconventional  Machines . ...  5 

Holl and/ SOLOMON. ..1 

ILLIAC  IV 3 

ILLIAC  IV/B65OO...I 

A  short  history  of  conventional  machine  design  will  preface 
a  discussion  of  the  antecedents  of  ILLIAC  IV,  namely  the 
Holland  and  SOLOMON  machines.   The  discussion  of  the  ILLIAC  IV 
organization  will  include  discussion  of  the  architecture  of 
the  B65OO. 


-16- 


II.  Operating  System  and  l/o  .  .  (Marv  Graham)  ...  2  classes 

(Tom  Mason  ) 

The  operating  system  is  presented  as  the  "glue"  binding 
the  B65OO  and  the  IKLIAC  IV.  The  l/O  subsystem  will  be 
explored  for  timing  and  operational  relationships. 


Ill  (a).   Introduction  and  General 

Overview  of  Software  .  .  (Prof.  Northcote)  .  .  1  class 

III  (b).   Software 11  classes 

Assembler/simulator  .  .  (Dave  Grothe   )  •  .  3 

(Gary  Grossman) 

GLYPNIR (Duncan  Lawrie)  .  .  k 

TRAIIQUIL (Norma  Abel    )  .  .  h 

The  major  languages  of  the  IIiLIAC  IV  are  presented.   An 
effort  will  be  made  to  highlight  the  useful  aspects. 

rV.   Communicating  with  the  User  .  .  (Pete  Alsberg)  .  .   1  class 
A  discussion  of  data  storage  devices  and  data  displays. 

V.   Applications   8  classes 

Numeric (Dave  Mclntyre  ) .    .    .    .   k 

Radar  (Morris  Knapp  ) ....  2 

Linear  Programming  .  .  (Serena  Yardley) .  ...  2 

The  major  uses  of  these  applications  will  be  sketched 
in  terms  of  how  ILLIAC  IV  is  essential  to  these  areas. 
Techniques  dependent  on  \inique  characteristics  of 
ILLIAC  IV  will  be  emphasized. 


3'7«2  Training  Programs 

As  the  delivery  of  ILLIAC  IV  draws  closer,  it  will  be  necessary 
to  initiate  training  programs  to  educate  people  in  the  use  of  ILLIAC  IV. 
Various  types  of  programs  are  currently  under  consideration  to  service 
three  types  of  people. 

(a)  Personnel  working  directly  for  the  ILLIAC  IV  Project 
at  the  University  of  Illinois  (U  of  l) 

(b)  Personnel  at  the  U  of  I  who  do  not  work  directly  for 
the  ILLIAC  IV  Project 

(c)  Personnel  outside  the  U  of  I  who  do  not  work  directly 
for  the  ILLIAC  IV  Project  but  are  interested  in 
becoming  users  of  ILLIAC  IV. 
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3. 7*3  Doc-umentation 

Due  to  the  growing  list  of  people,  "both  internal  and  external 
to  the  U  of  I,  who  utilize  ILLIAC  IV  Reports  and  Documents  and  the 
User's  Group- -FOURUM,  a  master  file  has  been  established.   This  master 
file  -will  effectively  automate  and  simplify  mailing  procedures  and  will 
centralize  user  descriptive  data  into  one  source. 

The  Library  of  ILLIAC  IV  documents  is  currently  undergoing 
re-organization  to  allow  for  future  expansion  and  to  make  it  easier  to 
use.   ACM  (Association  for  Computing  Machinery)  category  codes  will  be 
coded  to  existing  and  future  documents. 

3.7.^  ILLIAC  IV  Textbook 

Work  has  started  this  quarter  on  "The  ILLIAC  IV  Computer--A 
Programmer's  Text". 

The  book  is  an  attempt  to  present,  under  one  cover,  the  infor- 
mation necessary  to  instruct  a  programmer  in  the  use  of  ILLIAC  IV.   To 
this  end,  the  book  has  been  partitioned  into  two  distinct  sections.   The 
first  section,  the  chapters,  are  in  conventional  textbook  format 
including  problems  to  be  solved  by  the  student  and  topics  for  discussion. 
To  facilitate  the  learning  process,  the  information  presented  in  the 
chapters  is  quite  general  in  nature--only  those  concepts  necessary  to 
give  the  reader  a  basic  understanding  of  ILLIAC  IV  are  presented. 

The  second  section,  the  appendices,  refine  the  information 
presented  in  the  chapters.   The  appendices  act  as  a  reference  manual, 
or  a  handbook.   Either  section  of  the  book  can  be  read  independently 
of  the  other. 

Chapter  I  is  a  review--a  selected  history  of  computing 
machines.   The  evolution  of  digital  computers  is  described  in  terms  of 
the  problems  that  had  to  be  solved.   The  growing,  changing  capabilities 
of  computers  are  discussed,  leading  finally  to  ILLIAC  IV. 

Chapter  II  is  a  broad  outline  of  the  ILLIAC  IV  hardware 
structure.   A  sample  problem  (LaPlace's  equation  for  temperature 
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distribution  in  two  dimensions)  is  used  to  illustrate  how  the  components 
of  the  hardware  structure  are  tied  together. 

Chapter  III  describes  the  major  programming  languages  avail- 
able to  the  ILLIAC  IV  user.   The  sample  problem  of  Chapter  II  is 
re-visited  to  explore  the  use  of  the  assembly  language  ASK  and  two 
high  level  languages,  GLYPNIR  and  TRANQUIL. 

Chapter  IV  presents  the  system  software  that  is  available  to 
the  user.   The  l/o  system  including  the  B65OO  computer,  ILLIAC  IV  disk, 
DFC,  BIOM,  lOS,  CDC  and  TMU  is  described.   The  operating  system  is 
covered  in  a  skeletal  fashion  and  a  typical  job  is  traced  as  it  flows 
through  the  system.   Software  capabilities  of  specific  hardware  com- 
ponents such  as  remote  terminals,  graphics,  IMP  and  the  ARPA  network, 
and  various  other  peripheral  equipments  are  also  covered. 

Chapter  V  discusses  applications  of  ILLIAC  IV  and  indicates 
the  types  of  problems  that  are  especially  amenable  to  solution  on  such 
a  parallel  processor. 

3.8  Graphics 

The  following  specifications  are  currently  being  proposed  as 
the  graphic  requirements  for  ILLIAC  IV: 

I.   CRT  Requirements 

A.  20U8  X  2048  resolution  with  character  and  vector  capabilities. 

B.  Input  channel  of  ^10-^  bits/sec.   (see  II. B.) 

C.  Need  at  least  one  CRT  of  TV  size  for  viewing,  with  keyboard  for 
limited  interaction.   "Light  pen"  interaction  and  intermediate 
hard  copy  with  fair  quality  and  limited  quantity  (~15  sec/page) 
are  highly  desirable  if  not  too  costly. 

II.   Photographic  System 

A.  16  mm  and  35  nim  microfilm  cameras  with  sufficient  registration 
for  movies. 

B.  Required  rate  for  microfilm  estimated  at  1  frame/sec  with  an 
estimated  5OOO  instruc/frame  or  5OOO  characters/frame  (this 
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requires  CRT  input  channel  of  >  10^  bits/sec).   Development 
time  of  less  than  30  minutes  desired. 

C.   Final  hard  copy  from  microfilm  should  he  of  high  quality  and 
produced  at  high  rate. 

Vendors  are  currently  heing  contacted  to  submit  bids  on 
equipment  which  will  meet  the  above  specifications. 
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